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independent and identically A^(0, a^) distributed errors. Suppose that the parameter 

^ ' of interest is 6' = where a is a specified vector. Define the parameter r = 

^ ', c^P — t where the vector c and the number t are specified and a and c are linearly 

1 independent. Also suppose that we have uncertain prior information that r = 0. 
I ■ Kabaila and Giri (2009c) present a new frequentist 1 — a confidence interval for 

lO ' 9 that utilizes this prior information. This interval has expected length that (a) is 

^ ■ relatively small when the prior information about r is correct and (b) has a maximum 

. value that is not too large. It coincides with the standard 1 — a confidence interval 

^ ' (obtained by fitting the full model to the data) when the data strongly contradicts 
O ■ 

■ the prior information. At first sight, the computation of this new confidence interval 



^ ! seems to be infeasible. However, by the use of the various computational devices 
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and practicable. 
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1. Introduction 

Consider the linear regression model 

Y^X(3 + s 

where y is a random n-vector of responses, X is a known n x p matrix with lin- 
early independent columns, (3 — . . . ,f3p) is an unknown parameter vector and 
e ~ A^(0, (j'^In) where cr^ is an unknown positive parameter. Suppose that the pa- 
rameter of interest is ^ = (3 where a is specified p- vector (a ^ 0). Define the 
parameter r = j3 — t where the vector c and the number t are specified and a 
and c are linearly independent. Also suppose that previous experience with similar 
data sets and/or expert opinion and scientific background suggest that r = 0. In 
other words, suppose that we have uncertain prior information that r = 0. Ex- 
amples include having uncertain prior information that (a) one of the regression 
coefficients takes a specified value and (b) the linear regression consists of two 
parallel straight line regressions. "Higher order" terms in a linear regression model 
are often strong candidates for terms that could plausibly be zero. For example, 
for factorial experiments it is commonly believed that three-factor and higher order 
interactions are negligible. Indeed, this type of belief is the basis for the design of 
fractional factorial experiments. Another example is that it is commonly believed 
that the highest order terms in a univariate or multivariate polynomial regression 
are likely to be negligible. Our aim is to find a frequentist 1 — a confidence interval 
(i.e. a confidence interval whose coverage probability has infimum 1 — a) for 9 that 
utilizes this uncertain prior information, based on an observation of Y . 

One may attempt to utilize the uncertain prior information as follows. We carry 
out a preliminary test of the null hypothesis r = 0, against the alternative hypothesis 
T 7^ 0. We then find the confidence interval for 9, with nominal coverage 1 — a, based 
on the assumption that the selected model had been given to us a priori. It might 
be hoped that this confidence interval will have good coverage properties and an 
expected length that (a) is relatively small when the prior information is correct 
and (b) is not too large when the prior information happens to be incorrect. This 
assumption is false and, as pointed out by Kabaila (1995, 1998, 2005, 2009), Giri 
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and Kabaila (2008), Kabaila and Giri (2009b) and Kabaila and Leeb (2006), it leads 
to a confidence interval whose minimum coverage is typically far below 1 — a. In 
other words, this confidence interval fails abysmally to utilize the uncertain prior 
information. 

We assess a 1 — a confidence interval for 6 using the ratio (expected length of this 
confidence interval) /(expected length of standard 1 — a confidence interval). The 
standard 1 — a confidence interval is obtained by fitting the full model to the data. 
We call this ratio the scaled expected length of this confidence interval. Kabaila and 
Giri (2009c) describe a new 1 — a confidence interval for 9 that utilizes the prior 
information. This interval has scaled expected length that (a) is substantially smaller 
than 1 when the prior information that r = is correct and (b) has a maximum 
value that is not too much larger than 1. It coincides with the standard 1 — a 
confidence interval when the data strongly contradicts the prior information. This 
interval also has the attractive property that it has endpoints that are continuous 
functions of the data. 

Let © and f denote the least squares estimator of 9 and r respectively. Define the 
correlation coefficient p — corr(©,f). Also define the parameter 7 = T/-\/var(f). 
Both the coverage probability and the scaled expected length of the new 1 — a 
confidence interval are even functions of 7. An example of the performance of this 
confidence interval is shown in Figure 2 for the case that l — a = 0.95, n — p=l and 
p — OA. The top panel of this figure is a plot of the coverage probability of the new 
0.95 confidence interval for ^ as a function of 7. This plot shows that this coverage 
probability is 0.95 throughout the parameter space. The bottom panel of Figure 2 
is a plot of the square of the scaled expected length of this confidence interval as a 
function of 7. When the prior information is correct (i.e. 7 = 0), we gain since the 
square of the scaled expected length is substantially smaller than 1. The maximum 
value of the square of the scaled expected length is not too large. The new 0.95 
confidence interval for 9 coincides with the standard 1 — a confidence interval when 
the data strongly contradicts the prior information. This is refiected in Figure 2 
by the fact that the square of the scaled expected length approaches 1 as 7 ^ 00. 
All computations presented in the paper were performed with programs written in 
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MATLAB, using the Optimization and Statistics toolboxes. 

In Section 2, we describe the constrained minimization problem that needs to be 
solved to find this new confidence interval. To arrive at this description, Kabaila and 
Giri (2009c) have already used the following simplification techniques: (a) invariance 
arguments that take account of the form of the uncertain prior information, (b) a 
simply-implemented constraint on the new confidence interval that guarantees that 
it will coincide with the standard 1 — a confidence interval when the data strongly 
contradicts the prior information and (c) simplified expressions for the coverage 
probability and the criterion to be minimized. Even so, the coverage probability 
constraint portion of this minimization problem involves a continuum of constraints. 
Thus, at first sight, the computation of this new confidence interval seems to be in- 
feasible. In Section 3 we describe how this continuum of constraints can be replaced 
by a finite number or appropriately-chosen constraints. Even though this makes the 
computation of the new confidence interval feasible, a significant number of compu- 
tational issues remain to be solved. The solution to these computational issues is 
described by Giri (2008) and presented in detail in the present paper. We compute 
the double integrals for the coverage probability, scaled expected length and the cri- 
terion to be minimized by first truncating these integrals. In Section 4, we present 
bounds on the resulting truncation errors. In Section 5 we present some practical 
advice on how to make these computations work. In Section 6, we present a nu- 
merical example that illustrates the successful computation of the new confidence 
interval. 

2. Constrained minimization problem to be solved 

Let P denote the least squares estimator of /3. Let denote i.e. the least 
squares estimator of 9. Also, let f denote c^P — t i.e. the least squares estimator 
of T. Define the matrix V to be the covariance matrix of (0, f) divided by a^. Let 

Vij denote the th element of V. The standard 1 — a confidence interval for 6 
(obtained by fitting the full model to the data) is / = [0 — tn-p,i-f^/v^a, + 
tn-p,i-!^-\/viiO'] , where the quantile t^^a is defined by P{T < tm,a) — a ior T tm 
and a^^{Y - Xpf{Y - Xfi) /{n - p). 

We use the notation [a ± b\ for the interval [a — b,a -\- b\ {b > 0). Define the 
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following confidence interval for 9 



J{b, s) = — ^/vll^' b 



( 



T 



± ^/VllO' S 



( 



T 



where the functions b and s are required to satisfy the following restriction. 

Restriction 1 6 : R ^ R is an odd function and s : [0, oo) [0, oo). 

The motivation for this restriction is provided by the invariance arguments presented 
in Appendix A of Kabaila and Giri (2009c). We also require that the functions b 
and s satisfy the following restriction. 

Restriction 2 b and s are continuous functions. 

This implies that the endpoints of the confidence interval J{b, s) are continuous 
functions of the data. Finally, we require the confidence interval J{b, s) to coincide 
with the standard 1 — a confidence interval / when the data strongly contradict 
the prior information. The statistic |f|/ {a■^yv22) provides some indication of how far 
away T/((Ty^^^) is from 0. We therefore require that the functions b and s satisfy 
the following restriction. 

Restriction 3 b{x) = for all |a;| > d and s{x) = t„_p_i_« for all x > d where d is 
a (sufficiently large) specified positive number. 

Define p = corr(e,f) = V12/ ^/v^{v^, 7 = r/ A/var(f) = r/((T^/t^) and W = 
a /a. Let m = n — p. Note that W has the same distribution as y^Q/m where 
Q ~ Xm- Let fw denote the probability density function of W. Note that fwiw) — 
2mw fmimw^) for all w > 0, where denotes the Xm probability density function. 

For given 6, s and p, the coverage probability P(9 E J{b, s)) is a function of 7. 
We denote this coverage probability by c{'y;b,s,p). Part of our evaluation of the 
confidence interval J(6, s) consists of comparing it with the standard 1 — a confidence 
interval / using the criterion (expected length of J(6, s))/ (expected length of /). We 
call this the scaled expected length of J{b, s). This is an even function of 7, for given 
s. We denote this function by e{j; s). 

Our aim is to find functions b and s that satisfy Restrictions 1-3 and such that 
(a) the infimum of 0(7; b, s, p) over 7 is 1 — a and (b) 




(1) 
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is minimized, where the weight function v has been chosen to be 

v{x) = Ax + nix) for all a; G M, (2) 

where A is a specified nonnegative number and 7i is the unit step function defined 
by l-Lix) = for X < and 7i(x) = 1 for a; > 0. This weight function has also been 
used by Farchione and Kabaila (2008) and Kabaila and Giri (2009a). The larger 
the value of A, the smaller the relative weight given to minimizing e(7; s) for 7 = 0, 
as opposed to minimizing e(7; s) for other values of 7. For appropriately chosen A, 
the weight function ([2]) leads to a 1 — a confidence interval for Q that has expected 
length that (a) is relatively small when r = and (b) has maximum value that is 
not too large. 

The following theorem provides computationally convenient expressions for the 
coverage probability and scaled expected length of J(6, s). 

Theorem 1. (Kabaila and Giri (2009c)). 

{a) DeGne i{h, w) = h{h/w) w — s{\h\/w)w and u{h, w) = h{h/w) w + s{\h\lw)w. 
Also define \E'(x, ?/; /i, f ) = P{x < Z < y) for Z ~ N{fi,v). Now define the 
functions k{h,w,'y,p) = {i{h,w),u{h,w)] p{h — j),! — p'^) and k'^ {h,w,'y, p) = 
\E' [—tn-p,i-^w, tn-p,i-qw; p{h — 7), 1 — p^) . The coverage probability of J{b, s) is 
equal to 

/•oo pd 

(! — «)+/ / (^k{wx,w,'y, p) — k\wx,w,^, p)) (j){wx — 'y) dxw fwiw) dw (3) 

Jo J-d 

where denotes the N{0, 1) probability density function. For given h, s and p, 
0(7; 6, s, p) is an even function of 7. 

(b) The scaled expected length of J{b, s) is denoted e{j; s) and is equal to 

2 roc pd 

^^t 1 ^ E{W) J J (^(1^1) ~^«-p.i-f) '^(^^~^)^^^Vi^(«^)c?w^- (4) 



The method used to compute E{W) is described in Appendix A. Substituting 
into ([1]), we obtain that ([1]) is equal to 

foo fd 



tn-p,l-^ E{W) 



poo pa 

/ / (s(a;)-t„_p,i_|)(A + 
Jo Jo 



(j){wx))dxw fw{w)dw 



This is proportional to 

I'OO I'd 

I I {s{x) -tn-p,!-^) {X + (j){wx))dxw'^ fw{w)dw 
Jo Jo 

(I'd \ POO I'd 

/ s{x) dx — dtn-p^i-^\ + / / (ys{x) — tn-p^i-^) 4>{wx) dxvo^ fw{w) dw, 

(5) 

since EiW"^) = 1. Therefore, our aim is to find functions b and s that satisfy 
Restrictions 1-3 and such that is minimized with respect to the functions b and 
s, subject to the constraint that > 1 — a for all 7 > 0. 

Unless we specify parametric forms for the functions b and s, the computation 
of these functions (to solve the constrained minimization problem) will certainly 
be infeasible. So, we specify the following parametric forms for these functions. 
We require 6 to be a continuous function and so it is necessary that 6(0) = 0. 
Suppose that Xi, . . . ,Xq satisfy = Xi < X2 < ■ ■ ■ < Xg = d. Obviously, b{xi) = 0, 
b{xq) = and s{xq) = tn-p^i-q- The function b is fully specified by the vector 
(6(^2), . . . , ^(xg-i)) as follows. Because b is assumed to be an odd function, we 
know that b{—Xi) = —b{xi) for i = 2, . . . ,q. We specify the value of b{x) for any 
X G [—d, d] by cubic spline interpolation for these given function values. We specify 
the function s by the vector (s(xi), . . . , s(xq_i)) as follows. The value of s{x) for any 
X G [0, d] is specified by cubic spline interpolation for these given function values. 
We call xi,X2, ■ ■ -Xq the knots. We have taken these knots to be equally spaced. 

To conclude, the new 1 — a confidence interval for 6 that utilizes the uncertain 
prior information that r = is obtained as follows. Theoretically, the performance of 
the new confidence interval will improve as d increases and the spacing between the 
knots Xi decreases. However, the computation of this confidence interval becomes 
numerically unstable if d is too large and/or the number of knots is too large. So, for 
each candidate value of the parameter A, we carry out the following computational 
procedure for judiciously-chosen sets of values of d and knots Xj. 

Computational Procedure 

Let z = (6(^2), . . . , 6(xq_i), s(xi), . . . , s(xg_i)) . Define f{z) to be the objective 
function ([5]), thought of as a function of z. Also define c(z;7) to be (1 — a)— (IHl), 
thought of as a function of z for given 7. Minimize f{z) subject to the constraints 
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that s{x) > for all x G [0, d] and the nonlinear coverage constraints that c{z; 7) < 
for all 7 > 0. Plot the coverage probability of J{b, s), as a function of 7 > 0. Also 
plot 6^(7; s), the square of the scaled expected length, as a function of 7 > 0. 

Based on these plots, and possibly on the strength of our prior information that 
T — 0, we choose appropriate values of A, d and knots Xi. Extensive guidelines for 
this choice are presented in Section 4 of Kabaila and Giri (2009c). The confidence 
interval corresponding to this choice is the new 1 — a confidence interval for 9. The 
focus of the present paper is how this Computational Procedure can be made feasible 
and practicable. 

3. Implementation of the coverage probability constraints 

At first sight, the continuum of nonlinear coverage constraints c{z; 7) < for 
all 7 > would seem to make the Computational Procedure infeasible. However, 
Restriction 3 implies that, for any reasonable choice of z, c{z;^) — > as 7 — > 00. 
Also, for any given value of z, c(z; 7) is a smooth function of 7 > 0. This suggests 
that this continuum of constraints can be replaced in the computations by the 
following finite set of constraints: c{z;j) < for every 7 G {0, A, 2A, . . . , M}, 
where A is a sufficiently small positive number and M is sufficiently large. It is easy 
to check numerically whether or not given values of A and M are adequate. If the 
graph of 0(7; 6, s, p) falls below 1 — a for some values of 7 > then this choice is 
inadequate. On the other hand, if 0(7; b, s, p) > 1 — a for all 7 > then this choice 
is adequate. This numerical check corresponds to the following easily-proved result. 

Lemma 1. Suppose that z* minimizes f{z) subject to the constraints that s{x) > 
for all X e [0, d] and the coverage constraints that c{z; 7) < for aJJ 7 > 0. Also, 
suppose that z' minimizes f{z) subject to the constraints that s{x) > for all x G 
[0, d] and the coverage constraints that c{z; 7) < for every 7 G {0, A, 2A, . . . , M}. 
Ifc{z';-f) < for aJJ 7 > then f{z') = f{z*). 

For the numerical example presented in Section 7, we chose A = 0.5 and M — 50. 
That this choice is adequate is clear from the plot of c(7;6, s, p), as a function of 
7 > 0, in the top panel of Figure 2. 
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The function s needs to satisfy the continuum of constraints s{x) > for all 
X G [0, d]. Similarly to the coverage probability constraints, these could be replaced 
by the following finite set of constraints: s{x) > for each x G {0,6,26, ... ,d} where 
(5 is a sufficiently small positive number. However, it was found that the constraints 
s{xi) > \tn-p,i-q for z = — 1 were not too restrictive and, in practice, 

guaranteed that s{x) > for all x G [0,d\. 

The constrained minimization problem is solved numerically using the MAT- 
LAB function fmincon. The starting value of z was chosen to correspond to 
the standard 1 — a confidence interval I. In other words, for this starting value, 
b{x2) = 0, . . . , = 0, s{xi) = tn-pA-f, s{xg-i) = t„-p,i-|. The "Medium- 

Scale Optimization" option for this function is used. This option uses a Sequential 
Quadratic Programming (SQP) method described in detail in the documentation 
for the Optimization toolbox. 

4. Bounds on the truncation errors 

The double integrals in ([3]), (jl]) and the second term on the right-hand-side of 
([5]) are evaluated as follows. These integrals are first truncated with respect to 
w, followed by numerical evaluation of the truncated double integrals using the 
MATLAB function dblquad. In this section, we derive bounds on the resulting 
truncation errors. We use c to denote the upper endpoint of the truncated integral 
with respect to w. 

Define the truncation error 

POD I'd 

ei = / / {k{wx,w,'^, p) — k\wx,w,'^, pj) 4){wx — dxw fw{w) dw 

Jo J-d 

c I'd 

I {k{wx, w, 7, p) — k\wx, w, 7, p)) (f){wx — 7) dxw fw{w) dw 
J-d 

00 I'd 

/ {k{wx, w, 7, p) — k'^{wx, w, 7, p)) (f){wx — 'y) dxw fwiw) dw. 
J-d 

As proved in Appendix B, |ei| < P{Q > mc^), where Q ~ Xm- 

To find bounds on the other truncation errors, we will use the following lemma, 
which is proved in Appendix C. 
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Lemma 2. 



J Wfw{w)dw = y— p^^m-^^ > ™^^)' ^^^^^ <5 ~ Xm+V 



Define the truncation error 

f oo pd 



poo pa 

62 = / / {s{\x\)-tn-p,l-^)(l){wX-'y)dxw'^fwiw)dw 
Jo J-d 

PC pd 

- / {si\x\)-tn-p,i-^)(l){wx--f)dxw'^ fwiw)dw 

Jo J-d 

J {s{\x\) - tn-p,i-^) (f){wx - ^) dx fw{w) dw. 
As proved in Appendix D, |e2| is bounded above by 

where Q ~ xL+i- 

Define the truncation error 

POO pd 

^3 = / / - in-p,i-f ) dxw'^ fw{w) dw 

Jo Jo 

{s{x) — tn-p^i-^) (l){wx) dxw^ fw{w) dw 



Jo 



{s{x) — tn-p,i-^) (p{wx) dxw^ fw{w) dw 

Ic Jo ' ^ 

As proved in Appendix E, |e3| is bounded above by 

max,>o \s{y) - A2" r(f + 1) 

2 V m ^XfT^ ^' 

where Q ~ Xm+i- 

We may, very conservatively, assume that maXj^>o \s{y) — tn-p,i-f | ^ 10 at or 
near the solution to the constrained minimization problem. For any given value 
of m, it is easy to compute the values of c such that these upper bounds on the 
magnitudes of the the truncation errors are equal to some small specified positive 
number. 
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5. Some practical advice 

It was found that the computation of the coverage probabihty ([2]), which entails 
the computation of 

/ / [k{wx,w,'y, p) — k^{wx,w,'y, p)) (j){wx — 'y) dxw fw{w) dw (6) 

Jo J-d 

using the MATLAB function dblquad, was inaccurate for small n — p. The reason 
for this was found numerically to be the following. For small n — p the integrand of 
([6]) is non-zero only for w very close to zero; elsewhere in the interval [0, c] it is very 
close to zero. As a result, dblquad may largely "miss" the non-zero values of the 
integrand, leading to the inaccurate computation of ([6]). Our pragmatic solution 
to this problem is as follows. If c > 3 then we perform two numerical integrations 
using dblquad. The first numerical integration evaluates 

/ / {k{wx,w,'y, p) — k\wx,w,'y, p)^ (j){wx — 'y) dxw fw{w) dw 

Jo J-d 

and the second numerical integration evaluates 

PC I'd 

I I {k{wx,w,'y,p) — k\wx,w,'y,p))(j){wx — 'y)dxwfw{w)dw 
J2 J~d 

These two evaluations are then added to obtain the computed value of ([6]). 

To help prevent the occasional instability in the computation of the solution to 
the constrained minimization problem, the following bounds were applied: —100 < 
b{xi) < 100 for i = 2, . . . , g — 1 and s{xi) < 200 for i = 1, . . . , g — 1. In some cases, it 
was found that the spline defining the function h had oscillations that were clearly 
spurious. These oscillations disappeared when the endpoint constraints b'{—d) = 
and b'{d) = were introduced. These endpoint constraints are now part of the 
computational method. 

The computation of the solution to the constrained minimization problem can 
be quite delicate, especially when n—p is small. In some cases the computation may 
not converge to the solution to the constrained minimization problem, as evidenced 
by poor coverage properties of the computed confidence interval and/or spurious 
oscillations in the values of b and s. In this case, the computed solution is used 
as the starting value for another computation of the solution to the constrained 
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minimization problem. This procedure usually leads to the successful computation 
of the solution to this minimization problem. 

6. Numerical example 

Kabaila and Giri (2009c) present an example of the new 1 — a confidence interval 
that utilizes the uncertain prior information, for the case that p — —l/y/2, n—p — 76 
and 1 — a = 0.95. In the present section we consider the more computationally 
challenging case that p = 0.4, n—p = 1 and 1 — a = 0.95. For each candidate value of 
the parameter A, we carried out the Computational Procedure (described in Section 
2) for judiciously-chosen sets of values of d and knots Xj. Using the guidelines for the 
choice of A, d and the knots Xi presented in Section 4 of Kabaila and Giri (2009c), we 
chose A = 0.2, d — 30 and the equidistant knots Xi at 0, (d/6), . . . , d. The resulting 
functions b and s, which specify the new 0.95 confidence interval for 9 that utilizes 
the uncertain prior information, are plotted in Figure 1. The performance of this 
confidence interval is shown in Figure 2. The top panel of this figure shows that 
the coverage probability of this confidence interval is 0.95 throughout the parameter 
space. The bottom panel of Figure 2 is a plot of the square of the scaled expected 
length of this confidence interval as a function of 7. When the prior information 
is correct (i.e. 7 = 0), we gain since the square of the scaled expected length 
is substantially smaller than 1. The maximum value of the square of the scaled 
expected length is not too large. The new 0.95 confidence interval for 9 coincides 
with the standard 1 — a confidence interval when the data strongly contradicts the 
prior information. This is refiected in Figure 2 by the fact that the square of the 
scaled expected length approaches 1 as 7 — > 00. 
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Figure 1: Plots of the functions b and s for the new confidence interval for 6 when 
p = 0.4, n — p = 1 and 1 — a = 0.95. These functions are obtained using d = 30, 
A = 0.2 and the equidistant knots Xj at 0,{d/Q), . . . ,d. 
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Figure 2: Plots of the coverage probability and 6^(7; s), the squared scaled expected 
length, as functions of 7 of the new 0.95 confidence interval for 6 when n — p = 1, 
p = 0.4. These plots are obtained using d = 30, X = 0.2 and the equidistant knots 
Xi at at 0, {d/6), . . . ,d. 
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Appendix A. Computation of E{W) 

Using the well-known formula for the moments of a random variable with a 
gamma distribution (see e.g. Casella and Berger (2002, p. 130)), it may be shown 
that 

^(p^) = ./l£(^±i)_ 

^ ^ V m r(f ) 

By 6.1.47 on p.257 of Abramowitz and Stegun (1965), E(W) ^ 1 as m ^ oo. 
However, when m is even moderately large, r(m/2) is extremely large. To avoid 
problems with overflow, wc flrst compute In (r(| + y)) and In (r(y)) by using the 
MATLAB function gammaln. We then find E{W) by computing 



exp(-|ln(f)+ln(r(| + f))-ln(r( 



Appendix B. Derivation of the bounds on the truncation error Ci 

By changing the variable of integration from x to h — wx in the inner integral 
that defines ei, we obtain 

POD fdw 

ei = / / {k{h,w,^,p) - k^{h,w,^,p))(f){h-^)dhfw{w)dw. 

J c J —dw 

It follows from the proof of Theorem 1 of Kabaila and Giri (2009c) that k{h, w, 7, p) 
and k^{h,w,j, p) are conditional probabilities and so they belong to [0, 1]. Hence 

fdw 

-1 < / ik(K w, 7, p) — kHh, w, 7, p) I mh — 'y)dh< 1. 



/aw 
{k{h, w, 7, p) - k\h, w, 7, p)) (j){h - 7) dh 
-dw 



Therefore 



/oo poo 
fw(w) dw <ei < J fw{w) dw. 

Remember, fwiw) — 2mwfm{rnw'^) for all w > 0, where fm denotes the proba- 
bility density function. Thus 



/oo /"OO 
fw{w) dw — I 2mw /ot(?7t.w^) dw. 



Changing the variable of integration from w to y = mw'^, we find that this integral 
is equal to 

poo 

/ fm{y)dy^P{Q >mc^), 
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where Q ~ x^. Hence |ei| < P{Q > mc^). 

Appendix C. Proof of Lemma 2 

Observe that 



poo 

}fw{w)dw —I w 2mw fmimw^) dw, 

J c 



W] ,, ^ , 

•J c 

where, as in Section 2, fm denotes the Xm probabihty density function. Changing 
the variable of integration to y — mw^, the right-hand-side becomes 

l_ ^(f + i) P I ((m+l)/2)-l) , 

m r(f) 7^,2 2(™+i)/2r(f + 1) 
The result follows from the fact that the integral in this expression is equal to 
P{Q >mc^), where Q ~ Xm+i- 

Appendix D. Derivation of the bounds on the truncation error 

By changing the variable of integration from x io h — wx in the inner integral 
that defines 62, we obtain 



62 = ^ y ^K^^ ~ ^n-p,i-f^ (t>{h - 7) dxwfw{w) dw. 



Thus, 



poo pdw 

|e2| < / / m.B.-K\s{y) -tn-p,i-^\(t>{h- ^)dxw fw{w)dw 

Jc J-dw 

roo pdw 

= max \s{y)- in-p,i-# I / / (t){h-^)dxw fw{w) dw 

Jc J-dw 

POO 

< max Is (7/) - t„_p,i_a| / wfw{w)dw. 
The result now follows from Lemma 2. 

Appendix E. Derivation of the bounds on the truncation error 63 

By changing the variable of integration from x to h — wx in the inner integral 
that defines 63, we obtain 

/oo pdw 
J \s{h/w) - tn-p,i-^\(l){h) dhw fw{w) dw 

POO pdw 

< / max\s{y) — tn-p (j){h) dhw fw{w) dw 
Jc Jq y>^ ' ^ 

maxj,>0 p „ , , , 

< 7\ / w}w[w)dw. 
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The result now follows from Lemma 2. 
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